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I. INTRODUCTION 



t-j- \ Since its foundation more than 30 years ago, one of the main goals of Quantum Chromodynamics [QCD] is to 
^| 1 understand the structure of hadrons via the interactions of its elementary degrees of freedom, the quarks and gluons. 
In order to approach this fundamental question, one of the most appropriate starting points is to consider Coulomb 
gauge. In this gauge, there is a natural connection to the would-be physical degrees of freedom [1], and an appealing 
picture of confinement exists: the so-called Gribov-Zwanziger scenario of confinement [1-3], whereby the temporal 
component of the gluon propagator provides for a long range confining force in the infrared, while the transversal spatial 
Oh 1 component is infrared suppressed (and therefore does not appear as an asymptotic state). In the last years, important 
progress has been made in investigating the Yang- Mills sector of the theory in Coulomb gauge: the variational method 
in the Hamiltonian approach [4-7], the Lagrange-based (Dyson-Schwinger) functional formalism [8, 9] and also lattice 
calculations [10, 11] (see also [12-15]). On the lattice, the temporal gluon propagator appears to be largely energy 
independent and behaves like 1/q 4 in the infrared, whereas the static spatial gluon propagator is found to be vanishing 
in the infrared, in agreement with the Gribov formula [10]. In the framework of the functional formalism it has been 
shown that the total color charge of the system is conserved and vanishing, and the well-known energy divergence 

j- problem in Coulomb gauge disappears [16]. Moreover, within this approach, the Slavnov- Taylor identities [17] have 

been derived and perturbative results have been provided [9, 18]. In the quark sector, perturbative results have been 
presented [19], and in the heavy mass limit the confining potential has been studied, both for two and three quark 
£^ \ systems [20, 21]. 

„h ■ In general, the underlying equation for the description of meson bound states is the two-body homogeneous Bethe- 
Salpeter equation. In the rainbow-ladder approximation, this equation has been successfully used to describe the 
properties of light mesons (see, for example [22, 23] and for a recent review [24]), where the driving mechanism 
is the chiral symmetry breaking. Beyond this approximation, models with dressed vertex contributions [25-30] 
and unquenching effects [31-33] have been considered, and more sophisticated numerical methods to solve both the 
homogeneous and the inhomogeneous Bethe-Salpeter equation have been recently developed [34]. In spite of this 
success, an exact derivation of the meson or diquark bound state energies via Green's functions techniques (i.e., 
using the Dyson-Schwinger equations as opposed to the Bethe-Salpeter equation) has not been yet reported. The 
difficulty stems from the fact that the (irreducible) interaction kernel contains higher order vertex functions which 
in general can not be calculated exactly. The connection between the singularities of the Green's functions and the 
bound states dates back to the original work of Gell-Mann and Low [35] , who gave a rigorous proof of the two-particle 
Bethe-Salpeter equation. It is, moreover, one of the longstanding problems in relativistic quantum field theory, and 
in particular phenomenological studies of meson and baryon states, to identify the physical solutions of the poles 
of the Green's functions and to separate them from possible unphysical ones. For an analysis of the appearance of 
unphysical poles in QED bound states, see for example Ref. [36] and references therein. 

Based on the heavy mass expansion underlying the Heavy Quark Effective Theory [HQET] [37, 38] and with the 
truncation of the Yang-Mills sector to include only dressed two-point functions, we study nonperturbatively quark- 
antiquark and diquark states using Green's functions techniques. By means of functional methods, we explicitly derive 
the Dyson-Schwinger equation for the quark 4-point Green's functions and give an exact, analytical solution. This 
will enable us to verify that bound states are related to the occurrence of physical poles in the Green's functions and 
hence we will be able to provide a direct connection between the homogeneous Bethe-Salpeter equation (considered 
previously in Ref. [20]) and the singularities of the Green's function, at least within the scheme considered here. 
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Moreover, within the employed approximation we will show that the physical and unphysical poles naturally separate. 

This paper is organized as follows. In Sec. II we review the functional formalism and derive the equations of motion, 
along with the relevant higher order functional derivatives. We present the explicit derivation of the Dyson-Schwinger 
equations for the 4-point Green's functions for quark- antiquark systems. In Sec. Ill we give a brief survey of the 
results obtained for heavy quark systems. After making an expansion of the generating functional of Coulomb gauge 
QCD in powers of 1/m, where m is the quark mass, and keeping only the leading term, we present the main steps 
in the derivation of the (anti)quark gap equation, and the corresponding temporal quark-gluon vertices. In Sec. IV 
we consider both the 1-particle irreducible and amputated connected 4-point quark- antiquark Green's functions in 
the heavy mass limit. In this approximation, analytical solutions for the Green's functions are provided and the 
singularity structure, in particular the separation of the physical and unphysical states, is discussed. Moreover, the 
correspondence with the bound state solution emerging from the homogeneous Bethe-Salpeter equation is emphasized. 
Sec. V is devoted to diquark systems. Again, the separation of the poles and the connection to the Bethe-Salpeter 
equation (including flavor structure) is discussed. A summary and some concluding remarks are given in Sec. VI. In 
the Appendix, the explicit derivation of the Slavnov- Taylor identity for the quark-2 gluon vertex is presented. 



II. DYSON-SCHWINGER EQUATIONS FOR 4-POINT FUNCTIONS 

We begin by briefly reviewing the functional formalism used to derive the basic equations that will lead to the Dyson- 
Schwinger equations and Slavnov- Taylor identities. Following the notations and conventions from [19, 20], we work 
in Minkowski space, with the metric — diag(l, —1). Roman letters denote spatial indices, greek letters 

(/j,, v . . . ) denote Lorentz indices and superscripts (a, 6, c . . . ) stand for color indices in the adjoint representation. 
Unless otherwise specified, the Dirac flavor, spinor and (fundamental) color indices are labelled with a common index 
(a, (3 . . . ). Also, configuration space coordinates may be denoted with subscript (x,y,z . . .) when no confusion arises. 
The Dirac 7 matrices satisfy {7^,7"} = 2g^ v . 

Let us consider the explicit quark contribution to the full QCD generating functional [19]: 

Z[x,x] = y cxp d A xq a (x) i-f°D +ij-D -m ^9/3 00 j 

xexpjz J d i x[x a (x)q a (x) +q a (x)x a (x)} + iSym^- (2.1) 

In the above, Sym represents the Yang-Mills part of the action and X>$ generically denotes the functional integration 
measure over all fields. q a denotes the full quark field, q a is the conjugate (or antiquark) field, and x QJ Xa are the 
corresponding sources. The temporal and spatial components of the covariant derivative (in the fundamental color 
representation) are given by 

D = d -igT a a a (x), 

D = V + igT a A a (x), (2.2) 

where A and a refer to the spatial and temporal components of the gluon field, respectively, and the source terms 
for A and a (with sources J and p, respectively) are implicitly included in Sym- The structure constants of the 
SU(N C ) group are denoted with f abc 1 and the Hermitian generators T a satisfy [T a ,T b ] = l j abc T c and are normalized 
via Tr(T°T b ) = S ab /2. Also note the Casimir factor associated with the quark gap equation: 

N 2 - 1 

C F = (2.3) 

The field equation of motion is derived from the generating functional, Eq. (2.1), and the observation that the 
integral of a total derivative vanishes, up to possible boundary terms. We assume that these terms do not contribute, 
although this is not obvious due to the Gribov problem [2] (see [8] and references therein), and obtain: 

/ V ® Siq \y) exp {^ J d4x7 la( x ) il°D + ij-D -m ^q p (x)+i J d*x [x a (x)q a (x) + q a (x)x a (x)] + iS Y m^ = 0. 

(2-4) 

So far, the generating functional Z[J] generates both connected and disconnected Green's functions. The generating 
functional of connected Green's functions is W[J], where W = In Z. We introduce a bracket notation for the functional 
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derivatives of W, such that for a generic source denoted by J a (in this case the index a refers to all the attributes of 
both the quark and gluon fields, including the type): 

=<tJ a > • (2.5) 

OlJ a 

The classical fields are (we use the same notation for the classical and quantum fields, as is standard): 

^) = ^/^Wexp^} = I^. (2 . 6) 

The generating functional of the proper (one particle irreducible) Green's function is the effective action T (function 
of the classical fields) and is defined via the Legendre transform of W[J]: 

= W[J] - iJ a <S> a . (2.7) 

In the above, we have used the common convention that the index a is either summed or integrated over as appropriate. 
This gives 

SW 5T 
<iJ a >= y— = <S> a and <i$ Q >= — — = -J a , (2.8) 

OlJ a 0l9 a 

where the same bracket notation for derivatives of T with respect to the classical fields is used - there is no confusion 
between the two sets of brackets since the derivatives with respect to the sources and fields are never mixed. 

We now present the quark field equation of motion in terms of proper Green's functions, arising from Eq. (2.4) (and 
from which the gap equation and the 4-point quark Green's functions will be derived): 

<iq a (x)> = -i i"/ a d Qx + ij ■ V x - m iqp{x) 

+ [gT c ^] ap [a c (x)q fj (x) + <zp c (x)zxp(x)>] [.^V]^ [AUx)qp(x)+ <zJ c k (x)z X ^x)>] . (2.9) 

We retain for the moment the spatial part in the above equation, although later on this will be truncated out within 
the employed approximation scheme. 

At this stage, it is useful to introduce some notation for multiple functional derivatives with respect to fields 
and sources, which will be later on used to derive the Dyson-Schwinger equations. Consider the following partial 
differentiation expressions arising from the Legendre transform, Eq. (2.7) (we omit the configuration space arguments 
for clarity): 

— — <X(J)> = -»%] <i9 p %9~x%J~X(J)>, (2.10a) 

__ <r($)> = is[y] <*j i8 *j 7 ><»$ 7 r($)>, (2.10b) 

Old p 

where S[-f] = ±1 accounts for the fact that the fields may be Grassmann-valued, i.e. a minus sign appears when the 
index 7 refers to the following contribution to the sum (over "type" 7) 

<. . . iqyXiXy . . .>, <. . . iXyXiq-y ■ ■ ■> ■ (2-11) 
Combining Eq. (2.10a) and Eq. (2.10b) recursively, one easily finds the following useful relationship: 

S 



<iJsiJ a >= —r]0sS[^f, k] <iJ ( 5«J K X« < I' K < I> ) g« < I' 7 ><2J 7 »J ct >, (2-12) 



where the factor r)ps — — 1 if fields/sources of the type 5, (3 anticommutc. Of course, this type of relationship is 
standard, but the notation allows one to follow the various multiple derivatives and keep track of the signs efficiently 
The quark gap equation is derived by taking the functional derivative of the quark field equation of motion (in 
configuration space), Eq. (2.9), with respect to iq 7 (z), and omitting the terms which will eventually vanish when the 
sources are set to zero: 



<iq a {x)iq 1 {z)> — 1 ij°d 0x + 17 ■ V x — m 5{x — z) 



a 7 



-J d*y5(x-y) [r^^y <^(^(-)> ^{V>X^)> 



(2.13) 
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where the (configuration space) tree-level quark-gluon vertices, obtained from the quark equation of motion Eq. (2.9), 
are given by 

rj£* = [^W (2.14) 

r£5w, = -[rfL,, (2.i5) 

(omitting the trivial 5- function configuration space dependence). We now use the formula Eq. (2.12) to calculate 
the functional derivatives appearing in the bracket and obtain (for simplicity, we write the spatial arguments of the 
external fields/sources as subscripts and omit those for the internal fields/sources that are implicitly summed and 
integrated over): 

<iPyiXpx> = S[X,k] <ipliJx><i<&\iq 1 zi§K><iJKiXi}x>, (2.16) 



Siq JZ 

-j— <iJ? y iXf)x> = S[\,k] <iJ* y iJx><i<i>\iq 7Z i<i> K ><iJ K iX0x> ■ (2.17) 

When sources are set to zero, the brackets on the right-hand side correspond to connected and proper Green's 
functions. We can then identify the field/source type "A" with the gluon (temporal for the upper expression, spatial 
for the lower) and "k" with the antiquark, and notice that in both cases S[X, n] = +1. Introducing our conventions 
and notations for the Fourier transform, we have for a general two-point function (connected or proper) which obeys 
translational invariance: 

<iJ av iJf) x > = j dkW aP {k)e- lk < y - x \ (2.18) 
<i$ ay i$ 0x > = [ dkT a p{k)e- %k <y- x \ (2.19) 



and for the three-point (proper vertex) function 

<i<S> ay i<S>p x i<S> lz >= J dh dk 2 dk 3 r aj87 (fci, fc 2 , k 3 )e- lkl -y-^- x - lk3 - z (27r) 4 S(k 1 + k 2 + k 3 ), (2.20) 

(similarly for higher n-point functions) where dk = d 4 k / (2tt) 4 . The propagator W a p(k) and proper (1PI) two-point 
function r a p(k) are related via the Legendre transform. For the quark propagator, we have the standard relation 

W qqai {k)T qqifl {k) =S aP . (2.21) 

We now insert the expressions Eqs. (2.16), (2.17) into Eq. (2.13). After Fourier transforming to momentum space 
and with the definitions Eqs. (2.18, 2.19, 2.20) we obtain the quark Dyson-Schwingcr (or gap) equation: 

iWfc) = rg^(fc) + / doj {r^ aj8 (fe, - u , u k)w lg p K (u)T\ qaivr (w, -k, k u)w? a {k u ) 

k)W- qq MT qqAK7j (uj, -k, k - Lj)W? Aij (k -w)}. (2.22) 

Let us now derive the Dyson-Schwinger equation for the one-particle irreducible (1PI) 4-point function. As an 
illustration of the functional differentiation techniques, we present the explicit derivation of the first term in this 
expression and notice that the rest of the terms follow from an identical calculation (as explained below, the A-tcrms 
will be eliminated within our truncation scheme and hence in order to keep things simple, in the following derivations 
we will drop them completely). The goal is to derive the 1PI four-point function <iq a iq^iq T iq n > , and hence on the 
right hand side we will need, for example, a term of the type 

diq vt biq TW biq lz y 

Functionally differentiating Eq. (2.16) with respect to iq TW and using the product rule, we obtain: 



5iq TW Siq^ 



<ip a yiXf>x> = -S[A,/c]| 



<iPylj\> 



Siq 

+ <zp^J A X«<I> A «q rto ig 7Z i$ K ><iJ K «X/3r> : > 
+VtjVtk <iPyiJ\Xi$\iq 7Z i$ K ,> 



<i<&\iq lz i§ K ><iJ K iXpx> 

<iJ K iX/3x> 



5iq T 



(2.23) 
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For clarity of presentation, we only retain the first term in the product (and denote the rest with dots). Again, we 
make use of the formula Eq. (2.12) and obtain: 



<tpyiXpx>= S[\,k,h,v] <ipyiJ fl ><i$ tl iq TW i$ v ><iJ v iJ\><i$\tq~ /z i$ K ,><iJ K ,ixpx> +■ 



8iq TW Siq lz 

A last functional derivative with respect to the quark field q v t gives 
S 3 



(2.24) 



5iq, qt diq TW 5iq 



<ip a yiXPx> = 5[A, k, /x, i/] 
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5iq vt 



<ipliJn> 



<i® ll iq TW i® v ><iJ v iJ\><i$\iq lz i® K ><iJ K iXpx> 



<tpyiJ tl ><t$^iq rl ttqTwi®u><iJ^J\><i®\iq-fzi®K><iJ^Xi3x> 

S 



5iq vt 



<lJ v lJ\> 



<i$\iq yz i$ K ><iJ K iXpx> 



+V V t <*PyiJn> <i^^q T wi^u> <iJ v iJ\> <i$\iq v tiq 1 zi$K> <iJ^Xf)x> 
+Vr,rVT 7 Vr, K <ip%iJ IJ: ><i$ l tiq TW i$ l/ ><iJ v iJx><i$\iqj Z i$ K > <iJ^Xl3x> 
+ .... (2.25) 
Again, we take only the first term from the above sum and as before we use the formula Eq. (2.12). We obtain: 
S 3 



5iq vt Siq TW 5tq. 
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<iPy*Xf3x> = -S[\,K,lX,V,E,S\<ip C l J lJ e ><l$ e iq r]t l$5><lJ5lJ t j l > 

x <i$piq TW i$„><iJ„iJx><i$\tq- tz i$ K ><iJ K ,iXi3x> +■■■■ 



(2.26) 



As for the gap equation, when sources/fields are set to zero, we can identify the possible contributing field types 
denoted by internal indices and determine S[X, re, p, v, e, 5). Straightening out the ordering of the external quark lines, 
the full Dyson-Schwinger equation for the proper (1PI) four-quark Green's function in configuration space reads 



<i-Q a x l< liziq TW iq v t> = [gi°T a ] a i3 J dy <5(a 



V) 



x { \<iX[i x iXK><iq K nizio\><ip c x ipf>\ [<iq TW tq v ia^><ix^Xe><iq e iq ri ti<Ts>< l P b » l Py>] 

- [<*Xi3x l xs><iqsiqr l ti< 7 s>< t P c s l pt>] [<iqTwiqv l<J i><iXvix\>^q^qyzivi><iplipl>] 

- [<i~Xf) x iXK> <iq^qizi<J L x><tp c x ipi> <q TW iqntiaii(Tl> <ip\ip a y >] 

+ [<ix / 3x t xs><iqsiqvtiv c s><ip L s' l pi>] [< i qTw i qizi<jii<ri><ip b x ipl>] 

+ [^X^X^^q^ziqrw^Vt^^PxVy^ 

- [<iXj3xiX K > <iq K tq 7Z iqxiq v t> <tq TW tq v ia b ^><ix^XX> <ip b ^Py>] 

- [<iXpx l Xn> <iq^q7ziq TW ' i q>.><ix\ixs> <iq s iq n ^o- b e > <ip b e ip a y >] 

+ [<iXi3 X ix K ><iqKiq"/zi(T c x><ip c \ipi>] [<iq TW iq^iaZ><ix^xs><tq8iq n ^ b e><' l P b e l Pl>] 

- \<^Xfix l X&><^q&nnt^l><iplipt>\ [^q^iq^a^^x^X^^q^zi^^P^Pp] } 

+ ... (2.27) 

where the dots represent the A vertex terms which are not considered here and we have already replaced the tree-level 
temporal quark-gluon vertex with its expression Eq. (2.14). The expression Eq. (2.27) is diagrammatically represented 
in Fig. 1 and the terms have been reordered such that the first term corresponds to the diagram (a) of Fig. 1, the 
second term corresponds to diagram (b) and so on. The term presented explicitly above corresponds to diagram (i). 
Also, notice the minus sign in the term corresponding to diagram (h) — this is the origin of the minus sign arising in 
the kernel of the homogeneous Bcthc-Salpeter equation considered in Ref. [20] . 

The next step is to derive the connection between the amputated connected 4-point quark-antiquark Green's function 
(used conventionally in the Bethe-Salpeter equation), and the 1PI Green's function. Starting with the identity 



<iq0iq 7 XiX-<iX a >-- 



(2.28) 
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(f) (g) ( h ) (i) (j) 

FIG. 1: Diagrammatic representation of the Dyson-Schwinger equation for the 1PI 4-point quark-antiquark Green's function. 
Blobs represent dressed proper (1PI) 4-point vertex, solid lines represent the quark propagator, springs denote either spatial 
(A) or temporal (a) gluon propagator and cross denotes the tree level quark-gluon vertex. Internal propagators and 1PI vertices 
are fully dressed. 




FIG. 2: Relation between the 1PI (dark blob) and amputated (shaded blob) 4-point Greens function for the quark-antiquark 
system. Internal propagators and 1PI vertices are fully dressed. 

we take two subsequent functional derivatives with respect to generic sources, 1J5 and iJ\. Using Eq. (2.10b) recur- 
sively, choosing 5, A as needed, we eventually obtain the standard expression (in terms of our notation and conventions) 
for the 4-point quark-antiquark connected Green's function, written in terms of 1PI Green's functions (and omitting 
the A- vertex contributions, which will not play a role here): 

<«X q «X<5*Xa z X»j> = <*Xa*X7><*XA i Xe><*9 7 i 9K*9e i 9/3>< i X K *X<5><*X/3*X t) > 

+ [<^x a l Xi><^^Q^%><^x^Xv>] [<iX\^XT><n T n^l><vc^xs>\ < i P a K ipl> 

- [<iXa l Xr><iq T i%ioi><ix^X5>] [<ix\iXe><iq £ m i °l><iX{}iXri>\ <%p a K ip h v > ■ (2.29) 

The relation Eq. (2.29) can be further simplified by introducing the fully amputated Green's function, i.e. dividing 
by the quark propagators (cut the quark legs), as shown in Fig. 2. Before we consider the explicit form of Eqs. (2.27, 
2.29), it is necessary to briefly recall the heavy mass expansion and the truncation scheme that will be applied in this 
study. This is the topic of the next section. 

III. HEAVY QUARK MASS EXPANSION 

Consider again the quark contribution to the full QCD generating functional, Eq. (2.1). In the following, we briefly 
describe the derivation of the quark (and antiquark) propagator, in the heavy mass limit and under the truncation of 
the Yang-Mills sector to include only the dressed two-point functions. The presentation follows the original work of 
Rcf. [20]. 
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The full quark field is decomposed according to the heavy quark transformation 

q a {x) = e- tmx °[h(x)+H{x)] a , K{x) = > [P+q(x)] a , H a (x) = e m *° [P-q(x)] a (3.1) 

(similarly for the antiquark field), where the two components h and H are introduced with the help of the spinor 
projectors 

P ± = l(l±7°). (3.2) 

This corresponds to a particular case of the heavy quark transform underlying HQET [37], but in the functional 
approach considered here it is simply an arbitrary decomposition that will turn out to be very useful in Coulomb 
gauge. As will shortly be discussed, this will lead to the suppression of the spatial gluon propagator and spatial 
quark-gluon vertex at leading order in the mass expansion. 

The quark fields, decomposed according to Eq. (3.1), are now inserted into the generating functional Eq. (2.1). 
Further, the //-fields are integrated out, and an expansion in the heavy quark mass is performed (throughout this 
work, we will follow the usual terminology and use the expression "mass expansion", instead of "expansion in the 
inverse mass"). At leading order, the generating functional reads: 

Z[x,x] = y"z>*expj t J d 4 xh a (x)[id Ox +gT a a a (x)] a0 hp(x) 

xexpj* J d 4 x[e- mx °x a (x)h a (x) + e mxo h a (x) X a(x)] + iSym} + O (1/m) , (3.3) 

where the temporal component of the covariant derivative D has been written explicitly. Notice that, as a consequence 
of the decomposition Eq. (3.1), the spin degrees of freedom have decoupled from the system. Importantly, in the above 
we have kept the full quark and antiquark sources (rather than the ones corresponding to the large components of 
the quark field h, as is usually done in HQET). This means that the methods of functional formalism are preserved 
at leading order in the mass expansion, and in turn this implies that we are allowed to use the full Dyson-Schwinger 
equations of QCD, i.e. the gap and 4-point Dyson-Schwinger equations of the last section, replacing, however, the 
kernels, propagators and vertices by their leading order in the mass expansion. 

The gap equation Eq. (2.22) is solved together with the Slavnov- Taylor identity. This is derived from the invariance 
of the QCD action under a time-dependent Gauss-BRST transform [20] , and in Coulomb gauge it has the following 
expression: 

fc3°r^(fcl,fc2,fc 3 ) = l^ qqAa0i (k 1 ,k 2 ,k 3 )T-J(-ks) 

fc 3 



+Tg qaS {k 1 ) f^fci + q , h - q Q ; k 2 ) + igT d 



S/3 



+ 



T d q - c Ak2 + qo,h-q :k l )-igT d Y qqSfj (-k 2 ) (3.4) 



where k\ + k 2 + k^, = 0, go is an arbitrary energy injection scale (arising from the noncovariance of Coulomb gauge [17]), 
T- c is the ghost proper two-point function, T q .^ cq and T q .^ cq are ghost-quark kernels associated with the Gauss-BRST 
transform. 

Since in the generating functional Eq. (3.3) the tree- level spatial quark gluon vertex appears at 0(1 /m) in the 
mass expansion (which is here neglected) , and furthermore, under our truncation the pure Yang-Mills vertices are also 
neglected, it follows that the Dyson-Schwinger equation for the nonperturbative spatial quark-gluon vertex provides 
the result that Tq q A ~ 0(1 /m) (see [20] for a complete discussion and justification of this truncation). Moreover, the 
ghost-quark kernels involve pure Yang-Mills vertices and can be also neglected. Thus, under our truncation scheme, 
the Slavnov- Taylor identity reduces to 

k 3 Tl qaa p(k 1 ,k 2 ,h)=r qqa s(k 1 )[igT d ] 5 p- [ lg T d } aS r qqSfj (-k 2 ) + 0(l/m). (3.5) 
We then insert this into Eq. (2.22), together with the tree- level quark proper two-point function 

r £/3 (*0 = *Safi [ko - m] + O (1/m) (3.6) 

and the tree level quark gluon vertex 

rS^fci, to, ks) - [ 9 T a ] a[i + O (1/m) (3.7) 
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that follow from the generating functional Eq. (3.3). The nonperturbative temporal gluon propagator reads [9]: 

W£(k) = 5 ab W a(T {k) = 8 ab i-D aa {k 2 ). (3.8) 



2 



k 

Lattice results [11] motivate that the dressing function D aa is largely independent of energy and moreover, D aa is 
infrared divergent and likely to behave as 1/k 2 for vanishing k 2 (we will only use the explicit form of D aa in the last 
step of the calculation). Putting all this together, we find the following solution to Eq. (2.22) for the heavy quark 
propagator: 

"^'*»- h-^t +M i +0(1/wl - (3 - 9) 

with the constant [du = d 3 uj/(2Tr) 3 ] 

Ir = ^ CF f^^l + oil,m). (3.10) 

The (implicit) regularization of I r is signaled by the symbol "r". When solving Eq. (2.22), the ordering of the 
integration is set such that the temporal integral is performed first, under the condition that the spatial integral is 
regularized and finite. With the solution Eq. (3.9), we return to the Slavnov- Taylor identity, Eq. (3.5), and find that 
nonpcrturbatively the temporal quark-gluon vertex remains bare: 

Tl qaap (h,k 2 , k 3 ) = [gT a ] ap + O (1/m) . (3.11) 

Let us now briefly review the properties of the propagator Eq. (3.9) (see also [20]). As opposed to the conventional 
QCD quark propagator, which possesses a pair of simple poles, in this case we only have a single pole in the complex 
fco-plane, and the explicit Feynman prescription becomes important. It then follows that the closed quark loops 
(constructed from virtual quark- ant iquark pairs) vanish due to the energy integration, and this implies that the 
theory is quenched in the heavy mass limit: 



/ 



dko = 0. (3.12) 



[fco — m — X r + is] [ho + po — m — T r + is] 



Notice also that the propagator Eq. (3.9) is diagonal in the outer product of the fundamental color, flavor and spinor 
spaces. This corresponds to the decoupling of the spin from the heavy quark system, due to the decomposition 
Eq. (3.1). Indeed, differs from the heavy quark tree-level propagator [37] only through the mass term, and 

this is due to the fact that we retain the sources of the full quark fields, while in HQET one uses the sources for 
the large /i-fields directly. Finally, we also stress that the position of the pole has no physical meaning (the quark 
can never be on-shell). Since the poles in the quark propagator are situated at infinity (this becomes apparent after 
the regularization in I r is removed assuming, as indicated by the lattice [11], that the temporal gluon propagator is 
infrared enhanced), one requires infinite energy to create a single quark from the vacuum or, if a hadronic system is 
considered, only the relative energy (derived from the Bethe-Salpeter equation) is important. Indeed, it is precisely 
the exact cancellation of the divergent constants that lead to the separation of physical and unphysical poles. 

As discussed in Ref. [20] (and visible in the single pole structure of the quark propagator), the mass expansion leads 
to the separation of the quark and antiquark states. Whilst virtual gg-loops cancel due to the energy integration, 
Eq. (3.12), physical bound states for qq pairs should still exist. This means that we have to consider the gap equation 
for the antiquark propagator (and its Feynman prescription) separately. The gap equation for the antiquark propagator 
is derived from the full QCD gap equation Eq. (2.22), by reversing the ordering of the quark and antiquark functional 
derivatives that lead to the quark Green's functions: 

-IW-fc) = -T q % a (-k)-J dco {T\- q ^{-k,w,k-u)W q - ql p{-u)T^ 

+r b qqAS7j (-k,cj, k - u)W q - q7 p(-u)T q % ai (-cj, k,u>~ k)Wf Aij {k - w)} . 

(3.13) 

Again, we apply our truncation scheme and the above equation reduces to: 

iWa(-fc) = r£L(-*) + / awr^ 47 (-fc,w,fc-a;)w^^ (3.14) 
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FIG. 3: Crossed ladder diagram that contributes to the 1PI 4-point Green's function. The upper line denotes the quark 
propagator, the lower one, the antiquark propagator and springs denote the temporal gluon propagator. 



In similar manner, we derive the Slavnov- Taylor identity for the antiquark-gluon vertex: 

-k° 3 T d q - q ^ a (k 2 ,k u k 3 ) = +T qWS (k 2 ) [ ig T d ] T Sa - [igT^T^i-k^ + Oil/m). (3.15) 
The following antiquark propagator is obtained, as a solution to Eq. (3.14): 

W qqa0 (k) = + O (1/m) . (3.16) 

[ko + m — ± r + l£\ 

The corresponding temporal antiquark-gluon vertex is given by: 



r a qq<7a0 (h,k2, k 3 ) = - [gT\ a + O (1/m) . (3.17) 

In the above, notice the Feynman prescription of the propagator, as well as the sign of the loop correction. This will 
have the consequence that the Bethe-Salpeter equation for the quark-antiquark states can be interpreted as a bound 
state equation. There, the quark and the antiquark do not create a virtual quark-antiquark pair (as in the case of 
the closed quark loops), but a system composed of two separate unphysical particles (i.e., they are not connected 
by a single primitive vertex). Moreover, in the Bethe-Salpeter equation the so-called crossed box contributions (i.e., 
nonplanar diagrams that contain any combinations of nontrivial interactions allowed within our truncation scheme) 
cancel - in fact, this type of cancellation will also appear in the 4-point Green's function. This is due to the temporal 
integration performed over multiple propagators with the same relative sign for the Feynman prescription (similar to 
Eq. (3.12), but in this case the terms originate from internal quark or antiquark propagators). As a consequence, the 
Bethe-Salpeter kernel reduces to the ladder truncation [20]. 



IV. SOLUTION IN THE HEAVY MASS LIMIT 



Having reviewed the properties of the various vertices and propagators under truncation and in the heavy mass 
limit, we can now return to the formula for the 4-point Green's function Eq. (2.27) (combined with Eq. (2.29)), and 
the corresponding diagrammatic representations Fig. 1, Fig. 2, and apply our truncation scheme in the heavy quark 
limit, at leading order in the mass expansion. 



A. One particle irreducible Green's function 



For the quark-antiquark system, we consider the flavor non-singlet Green's function in the s-channel, where the 
quark and the antiquark are regarded as two distinct flavors (but with equal masses). Hence, the diagrams (a), (c) 
and (i) of Fig. 1 are excluded. The diagram (b) (crossed ladder type exchange diagram) explicitly reads in momentum 
space and without the minus sign prefactor (see also Fig. 3): 



/ 



-,(0)a 
qqaaS 



ipi, -Pi - u,u))W qqSlt> (pi +uj)T b q + (pi +w,p 4 ,-pi -P4-uj 



m g ^(P3, W -P3, -w)WWa(P3 - w)r^A 7 (P3 ~ W,P2,Pl + Pi + w)] W^(-Cj)W^( Pl + p A + u) (4.1) 



10 




- - CD 




»+p 1 +p 4 - 




FIG. 4: Truncated Dyson-Schwinger equation for the 1PI 4-point Green's function in the s-channel. Same conventions as in 
Fig. 1 apply. 



The above expression refers to the full QCD case (but omitting the spatial Waa contributions). The lower line 
corresponds to the antiquark, so we first rewrite it in terms of the (full QCD) antiquark propagators and vertices, 
using W qqtlX (P3 ~ <*>) = [- w qq(u ~ P3)] ^ and ^qqar^ = [- T gq<r\ Tfi - Tncn we a PP!y the mass expansion by inserting 
the appropriate propagators Eqs. (3.9, 3.16) and vertices Eqs. (3.11, 3.17). Noticing the energy independence of the 
temporal gluon propagator, one sees then immediately that the energy integral vanishes, due to the fact that both 
quark and antiquark propagators have the same Feynman prescription relative to uj : 



J dLJ Wq q (pi + U})W qi} (u} - p 3 ) - / du 



1 



W 



p\ 



■ ie 



U} 



Pi 



■ ie 



= 0. 



(4.2) 



That this crossed ladder diagram vanishes within the scheme considered here is exactly the analogous to the absence 
of the crossed ladder contribution to the kernel of the homogeneous Bethe-Salpeter equation in Ref. [20]. 

Turning to diagram (d), this contribution involves the Tq qaa vertex. In the Appendix it is shown that under 
the scheme considered here, the Slavnov- Taylor identity furnishes the result that the nonperturbative quark-2 gluon 
vertex (and hence the diagram (d)) vanishes. Again, this result is analogous to the discussion of the kernel in the 
homogeneous Bethe-Salpeter equation from [20]. 

Let us for the moment discard the diagrams (f) and (g), which include the 1PI 4-point quark Green's function, and 
the diagram (e), containing a 4 quark- gluon vertex. Then we are only left with the diagram (h) and the rainbow-ladder 
term (j). This simplification enables us to derive a solution for the corresponding (truncated) equation for the 1PI 
4-point quark Green's function. With this result at hand, we will then return to the diagrams (f), (g) and (e), and 
explicitly show that they cancel (and hence our assumption is justified). With these observations, Eq. (2.27) reduces 
to: 



- [9T a l°] afj j dy5(x - y) { [<ix^ x %X&><n5nr ) ^o%><ip%ip d v >\ 



(4.3) 



+ ^Xj3x l Xu><iq,M^Tw i qvt>< l x^XK><iq K iq lz i^ b x><ip h x ip a y>} . 

It is convenient to express the resulting equation in momentum space. We define 
<iq ax iq lz iq TW iq v t>= J dk 2 dk 3 dk 4 e^-x^k^k^k^i)^^ ^ ^ ^^(k, + ^ + ^ + ^ (4 4) 

and arrive at the following Dyson-Schwinger equation for the 1PI 4-point quark Green's function in the s-channel 
(shown diagrammatically in Fig. 4): 



ri 4 2 T „(pi,p 2 ,P3,P4) = - / duj 



r (0)o , 
1 qqaa/3 VP 1 ' 



-Pi - UJ,Uj)Wqqf38{pi + Uj)r C qq aSn { Pl + UJ,p 4 , ~ Pl - p 4 - U>) 



X [ T iqar t ,(P3,P2 ~ U,Pl + Pi + Oj) Wq qftK (u - p 2 )T\ qaK1 (u - P 2, Pi , -U>)] Wf a (-w) W% (pi + Pi + U}) 



I 



duT { ^„, a ( Pl , -pi - U,Uj)Wq ql3u (p 1 +W)r^] r?) (pi + W,p 2 - U>, f> 3 , Pi) W qq ^ K ((J - p 2 )V\ q(Tt . 1 (w - P 2 , Pi , ~U)) 



qqcrafi * 



(4.5) 



In order to proceed, we make the following assumption for the energy and momentum dependence of the function 

r (4) : 



r( 4 )( Pl ,p 2 ,p 3 ,P4) = r( 4 )(p ;pi 



(4.6) 
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FIG. 5: Momentum routing for the diagram (g). See text for details. 



with Pq = pi + p 2 - This implies that in the above equation the 4-point function T$, TV (pi + u,p 2 — cj,p 3 ,p 4 ) does 
not depend on the integration variable w , and hence we can separate the energy and three-momentum integrals. 
Moreover, with this Ansatz we will be allowed to Fourier transform the resulting spatial integral back to coordinate 
space, as shall be explained shortly below. 

We identify the antiquark component of Eq. (4.5) (lower line of Fig. 4) and insert the expressions Eqs. (3.9, 3.16), 
for the quark and antiquark propagators, along with the vertices Eqs. (3.11, 3.17). We also note the Fierz identity 
for the generators 



2 [T a U t T % = S - 



j-j-S a0 Ssj. 



(4.7) 



After completing the energy integration and with the definition Eq. (3.8) for the temporal gluon propagator, Eq. (4.5) 
simplifies to : 



[P - 22 r + 2ie] r£) (P ;pi + Pi) = 



9 



+i- 



i 

N. 



2 varjOr-y 



j dcjW aa (io)W aa (pi +p4 + uj) 
J dQW aa (uS)r%> TV (Po ; pi + Pi + to) . 



Let us now make the following color decomposition for the function r( 4 ): 



r(4) 



X X _|_ A A 



(4) 



(4.8) 



(4.9) 



where and T 2 4 ^ are scalar functions. At this point, it is convenient to Fourier transform back to coordinate space. 
In general, since Eq. (4.5) might in principle contain momentum-dependent vertex functions, as well as mixing of 
energy and three-momentum variables, this transformation could not be carried out. However, in our case (Eq. (4.8)) 
momentum-dependent vertices are absent and moreover, with the Ansatz Eq. (4.6), the energy and three-momentum 
integrals have separated such that the spatial integral is performed only over two functions (spatial gluon propagator 
and the spatial component of the quark 4-point function, both functions of momentum squared) . Hence the Fourier 
transform simplifies to: 



/ 



ctQW aa {Q 2 )Y^{P f) -{q + uj) 2 )= / dxe-^ 3 W aa {x)T^{P ;x), 

with x = \x\. Sorting out the color factors, it is straightforward to obtain for the components T 2 4 ^: 
P (4), p , _ ( 9 2 \ 2 W aa (xfN c [(P - 2X r ){Nl -2) + ig 2 C F W^{x)] 

-(4) 



(4.10) 



T^(P ;x) = i 



2NJ [ Po _ 21 r + i^-W^x) + 2ie][P Q - 21 r - ig^C F W a<J {x) + 2ie] 

2 \ 2 Uf („\2 



2N C 



w aa { x y 



P -21 r + i^ w -W aa {x)+2ie 



(4.11) 



where x is the separation associated with the momentum p\ + j>4- Inserting the above results into the decomposition 
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FIG. 6: Lowest order in the perturbative expansion of the proper four quark-gluon vertex. See text for details. 
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FIG. 7: Perturbative expansion of the diagram (e). 
See text for details. 



Boxes comprise m and n internal gluon lines, respectively, with m,n> 1. 



Eq. (4.9), we find the final formula for the 1PI quark Green's function: 

g 2 V W aa (xf 



^a-yrri [Po'l X ) 



2iV c J P - 21 r + i^W aa (x) + lie 



(Po - 21 r )N c (N 2 c - 2) + ig 2 N c C F W aa (x) 
P - 21 r - ig 2 C F W aa (x) + 2ie 



(4.12) 



Having derived the solution Eq. (4.12) for the 1PI Green's function, we return to the diagrams (f), (g) and (e) 
and show that they do not contribute to the final result. To see this, we first consider the diagram (g) and study 
the energy integral (see also Fig. 5). Noticing that the energy dependence of the internal four-point function can be 
written as 



r( 4 )(p + ^o) 



[oj + X + ie] n ' 

where X is a combination of constants, n = 1, 2 and m = 0, 1, the energy integral takes the form 



1 



[luq + Xi + le] 



= 0. 



(4.13) 



(4.14) 



Clearly, this integral is a generalization of Eq. (3.12) and this vanishes, just as for the loop corrections in the kernel of 
the Bethe-Salpeter equation from Ref. [20] and the diagram (b) from above. An identical calculation for the diagram 
(f), recalling that the lower line corresponds to an antiquark propagator, leads us to the fact that this integral is also 
vanishing. 

Finally, we are now in the position to show that the diagram (e), containing the 4 quark-gluon vertex, is also 
vanishing. The argumentation is based on our previous findings, namely that the diagrams (f) and (g), containing 
the 1P1 quark Green's functions, are zero. Because we exclude the Yang-Mills vertices, the external gluon of the 
five-point (proper) vertex can only couple to the (upper) quark or (lower) antiquark line. Then the 4 quark-gluon 
vertex at lowest order contains the diagrams shown in Fig. 6. The crossed diagrams vanish, based on integrals of the 
type Eq. (3.12) or Eq. (4.14), and in fact this can be generalized to the case with an arbitrary number of external 
gluon legs. The remaining contributions can then be included in the diagram (e) and can be written as a combination 
of diagrams of the form shown in Fig. 7, where the boxes contain an arbitrary number of internal gluon lines (ladder 
resummation). On the other hand, as a result of the Dyson-Schwinger equation for the 1PI quark Green's functions 
(as in Fig. 4), the 4-point functions contained in the diagrams (f) and (g) can also be written as a ladder resummation, 
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FIG. 8: Truncated Dyson- Schwinger equation for the fully amputated quark-antiquark 4-point Green's function in the s-channel. 
See text for details. 

and with the explicit gluon lines, diagrams (f), (g) coincide precisely with the internal substructure of the two terms 
in the diagram Fig. 7. Hence, the perturbative series of diagram (e) has been reorganized such that although the 
function T qqqqa itself does not vanish, this 5-point interaction vertex and the gluon line on top of it do indeed vanish 
at every order perturbatively. In turn, this implies that our original assumption is correct and the solution Eq. (4.12) 
is valid at every order in perturbation theory. As a side remark, we note that the nonvanishing of the 5-point function 
relates to the existence of three-quark bound states in the Faddeev equation, in the ladder approximation [21]. 



B. Amputated Green's function 



In the following, we consider the Dyson-Schwinger equation for the fully amputated 4-point quark-antiquark Green's 
function in the s-channel (which we denote by G^) and, with the simplifications outlined in the previous Section, we 
will derive a solution to this equation. The motivation for studying the fully amputated Green's function in addition 
to the proper function is the fact that the homogeneous Bethe-Salpeter equation is based on the amputated 4-point 
function. Since the reduction from Eq. (2.27) to Eq. (4.3) is valid for the 1PI Green's function, this must also hold 
for the amputated Green's function. We will explicitly verify that this is satisfied, and will analyze the position of 
the poles and show that they relate to the Bethe-Salpeter equation for physical states. 

The appropriate Dyson-Schwinger equation for the fully amputated 4-point quark-antiquark Green's function can 
thus be obtained from the formula Eq. (4.3), by replacing the 1PI Green's function T^ 4 ^ with the amputated Green's 
function G^> using the expression Eq. (2.29) and cutting the legs. This equation is the inhomogeneous ladder Bethe- 
Salpcter equation and writing it in terms of the usual variables reads, in momentum space (see also Fig. 8): 



,6 "I 

qqa\ T1 



Gi% TV (p + ,p-:k + ,k_) = w£(p-k)[T a qqry ] av [T\ 

-J dq \rfjw qq {q + )] ^ [W qq {q.)T b qq<J ]^ W%[p- $)G$. T „((f + , g_; k+, fc_).(4.15) 

In the above, the momenta of the quarks are given by p + = p + £P, p- = p — (1 — £)P (similarly for k and q), and 
£ is the momentum sharing fraction. P indicates the dependence on the total four momentum, which will become 
important for the investigation of the bound-state contributions to the Green's function. It comes as no surprise that 
the inhomogeneous Bethe-Salpeter equation is already in the ladder truncation, similar to Eq. (4.3). 

The right hand side of equation Eq. (4.15) does not depend on the external energy p , implying that the 4-point 
function G^ (p + , p- ; k + , fc_) has to be independent on the relative energies po,k , and we further assume that G^ 
depends on the relative spatial momentum p~k. We then identify the antiquark line as before, replace the heavy 
quark and antiquark propagators and vertices with the expressions Eqs. (3.9, 3.16, 3.11, 3.17), and perform the energy 
integration. We arrive at the following expression: 

Gi% TV (P ;p- k) = g 2 T« v T? 7 W aa (p- k) + g 2 T^ - ± - — J dqW^(p-q)G^ TV (P ;q-k). (4.16) 

Having integrated out the energy, it is convenient to rewrite the above formula back into coordinate space. Using 
the definition Eq. (4.4) and noting that the Fourier transform is as before given via Eq. (4.10), the equation Eq. (4.16) 
simplifies to (the separation x corresponds to the momentum p — k): 

GW rv (P ; x) = g 2 T a an T^W aa {x) + 9 2 T^ - ± - — W^(x)G^ Tri (P ;x). (4.17) 
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Again, we make a color decomposition of the function G^ : 



,(4) 



(4.18) 



(G^ and G^' are scalar functions), use the Fierz identity, Eq. (4.7), to sort out the color factors, and obtain the 
following results for the components G^\ G^: 



.(4) 



G?(P ;x) = [2- 



(P -21 r ) 2 W aa (x) 



2 J [P Q - 21 r + ifcW a<r (x) + 2ie][P - 22 r - i\ (N c - ij W aa {x) 
G?(P ;x) = -(£) Cftzj^gWg) . 

Replacing these results in the formula Eq. (4.18), we get the final result for the function G^: 



2ie\ 



(4.19) 
(4.20) 



Gi% Tri (Po;x) 



(P -21 r )W aa (x) 



2 P - 21 r 



2 \ , 



■W aa (x) + 2ie 



(Po - 2X r 



P - 2X r - ig 2 C F W <J<J (x) + 2ie 



(4.21) 



A direct calculation shows that our result for the amputated 4-point function is related to the result Eq. (4.12) for 
the 1PI Green's function, via the formula Eq. (2.29) in the s-channel (or alternatively, Fig. 2): 



it 4 ) r > + ,p-,fc_,fc + ) = G( 4 ) ;T > + , P -;fc + ,fc_) - wf a (v- k) [r| i [rg 



(4.22) 



A few comments regarding the structure of the above equation are in order. Notice that despite the truncation, the 
denominator structure of the 1PI and amputated Green's functions is identical in both color channels. Moreover, in 
this approach the physical and nonphysical poles disentangle automatically. Using the form Eq. (3.8) for the temporal 
gluon propagator, the denominator factor of the color singlet channel can be rewritten in the form 



1, 



dcoD aa (6j) 



C F [1 



(4.23) 



and hence the bound state (infrared confining) energy Po res (x) = a\x\ emerges as the finite pole position of the 
resonant component (first term in the bracket of Eq. (4.21) or Eq. (4.12)), for arbitrary number of colors, assuming 
that D aa (uj 2 ) ~ 1/uj 2 as indicated on the lattice [11]. This provides an explicit analytical dependence of the 4-point 
Green's function on the qq bound state energy, which results from the homogeneous Bethe-Salpeter equation from 
Ref. [20]. The: overall denominator factor for Eq. (4.21) is part of the normalization and has the explicit form 



Po 



du D aa (uj) 



C, 



2N r 



(4.24) 



which (like the quark propagator) implies a pole position at infinity when the regularization is removed. This factor 
does not appear in the homogeneous Bethe-Salpeter equation. Notice that the result, Eq. (4.21), is independent 
of £ (the momentum sharing fraction). That the physical pole position should be independent of £ is clear (c.f. 
Refs. [20, 23]); that the rest is independent of £ is presumably a consequence of the truncation scheme considered 
here. 



V. 4-POINT GREEN'S FUNCTIONS FOR DIQUARKS 



Let us now consider the diquark 4-point Green's function. This can be easily obtained from the equation Eq. (2.27) 
for quark-antiquark systems, by interchanging the quark legs and inserting the appropriate minus signs. We obtain 
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(a) 



(b) 



(c) 



(d) 



(e) 





(f) 



(g) 



(h) 



(i) 



(j) 



FIG. 9: Diagrammatic representation of the Dyson-Schwinger equation for the 1PI 4-point diquark Green's function. Same 
conventions as in Fig. 1 apply. 



(see also the diagrammatic representation from Fig. 9): 

<iq ax iq TW iq 7Z iqr,t>= [gj°T a } af 3 J dy5(x - y) 

x {- [<*X/3x*X«><^5 7 ^><zp^>] [<iq TW iq^iai><ix^X5><iqsn v ^ b e><iP b s iPy>] 

+ [<*X/3 a; *XK><^ K ^72 iCr A><*PA^><9ru.*'?r,t«Cr^Cr^><l/?| , 1 Zp^>] 

- \< % X^XK><nK^zi<J\><ip\ipi>\ \<iq TW iq v io\><ix^Xe><nen^i>^p\^p a y >\ 

- [<*x>x*Xk><*5kW?7z^><%™^^ 

+ [^X^XSXiqs^vt^tXV^P^] [<*9ru. i ^* cr ^><% i XA><«gA i 97^ iCr K ><2 P^Py > ] 

- [<*X/3 :c i Xft><«g K »g ru ,«'7 7z «'?A><«XA*X<5><% i 9')t* cr £>< i PeVy>] 



> 



+ [<*X/3o; i XK><^K^ru.*'?72 l ^t (7 A>< ;i PA*Py 

- [<ix~{)x % x&><nm^e%><wlwt>\ [<iq TW iq~{zi<Tiicrl><ip b x ipl>] 

+ [<iX0 X iX6><iq~6 i qr)ti(7 c s ><ip c e ipt>] [^q TW iq^ (j i><ix^XK><iq K iqizi(?\><ip\ipl>\ } + .... 



(5.1) 



As in the qq case, the dots represent the A vertex terms which are not considered here and the tree-level temporal 
quark-gluon vertex has been replaced with its expression Eq. (2.14). Also, notice that the diquark is antisymmetric 
under the exchange of two quark legs (a f> t or 7 H rj). We must therefore explicitly keep track of which quark is 
which and thus introduce flavor structure with an arbitrary number of equal mass flavors (before, we had to distinguish 
between quark and antiquark). 

As before, we analyze the diagrammatic representation from Fig. 9 and show that the same type of cancellations 
occur. Starting with the diagram (a) in momentum space, we notice that this is a crossed ladder type exchange 
diagram (see also Fig. 10): 



X 



± — ;roi — i 1 ro j 

[r|,„ M (P2,P4 + w,Pi +P3 -u)W 9gil s(-p i -u)T b 9gaSll (-p i -u,p i ,u)] W°£(w)W*(pi +p 3 -w). (5.2) 

It has been already shown that the integral over the quark propagators (with the same Feynman prescription) vanishes 
and thus the diagram (a) is zero. A similar type of integral arises in the diagram (j), and hence this term is also 
not giving a contribution. Further, the diagrams (b) and (h) are zero, since the corresponding quark-2 gluon vertex 
vanishes according to the Slavnov- Taylor identity, Eq. (A. 14). As in the case of the qq system, let us for the moment 
assume that the integrals (f), (i), containing the diquark 4-point vertex, and (g), containing the 5-point function 
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Pj+Pj-w 

FIG. 10: Crossed ladder diagram contributing to the 1PI 4-point Green's function in the diquark channel. 



Yqqqqa are also zero, and solve the Dyson-Schwinger equation with the remaining terms. Having derived the solution, 
we will then return to the diagrams (f),(i) and (g) and show that they vanish. The Dyson-Schwinger equation for 
the diquark 4-point function with the remaining diagrams, i.e. diagrams (c), (d) and (e) explicitly reads, in 

momentum space: 



^arjrj (Pi , P2,P3,Pi) 



r ^ia/?(Pl. ^ ~ Pl^)Wq q[3K {pi + Oj)T C qq(TK7 ( Pl +UJ, P3 ,-UJ-p 1 -p 3 ) 



X [r b q qaTV { P 2, W - P2, -u)Wqq UE (p 2 - w)!^^ ( P2 - UJ,p 4 , Pl + P3 + W )] W£{~Uj)W£{ P l + P 3 + w) 



+ 



/ 



T £q?aP (Pi . - Pi, W) Wqgffsipi + ^qqaSr, {Pi + U ,PA, ~U ~ Pi ~ Pi) 

x [^^(pa.w -p 2 , -w)W ftI ,A(P2 - Lo)T c qqaXl {p 2 - uj, P3 , Pi + P4 + u)] W£{-w)W^{pi + Pi + ui) 

- J duJ r ^t%(Pl. ~ w "Pi' u)Wq qf) x{pi + W) [r| 9ffTi ,(p 2 ,W -P2, -u)Wq qvK (p 2 - W)] 
X r it?r,(Pl +^,P2 - W,p 3 ,P4)W^(-w). 



(5.3) 



Inserting the expressions for the quark propagator Eq. (3.9) and vertex Eq. (3.11), resolving the color factors using 
the Fierz identity, Eq. (4.7), and simplifying the energy integrals, we have (noting the flavor channels with the §f 
factors): 



r ar-l,(Pl,P2,P3,P4) 



1 



P° - 2m - 21 r + 2ie 

.2 x 2 



9 

2N r 



[(7V C 2 + 1) 8 ai 8 Tri - 2N c S ari 8 T1 ] J dQ W aa {j3)W aa {pi +p 3 + Q) 

+ Z (^V~) 5 f ari 8l 1 [{Nl + l)8 av 8 T1 -2N c 8 ai 8 Tn ] j duj W aa {uj)W aa (pi + p 4 + uj) 
+ -^-$i\ s L [N c 8 aK S TX - S aX 8 TK ] J du 



ujq + pi — m — Zr + is ujq — p\ + m + X r — is 
x J dujW aa {Q)T^ n { Pl +Lo,p 2 -Lo,p 3 ,p 4 )Y (5.4) 

Noticing the distinctive energy, momentum and flavor structure of the equation, a suitable Ansatz for the diquark 
function is 

^ OLT^T] (Pi , P2,P3,Pi) ;pi +P3> + ;pi+P4) (5.5) 

where the superscripts c and p stand for the crossed and parallel configurations, respectively, Pq = p\ + p 2 as before, 
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and TaT^Ti still retain the color indices. With these notations, Eq. (5.4) decouples via the flavor structure to: 

7 2 1 



(p 0] p 1+ p 3 ) = (-i) 



2N C P° -2m- 2Xr + 2ie 

x 



| [(N* + 1) 6 aj S TV - 2N c S ar] S T7 ] J duW aa {ui)W aa (ft + p 3 + <3) 
+ [N c 5 aK S TX - S aX 5 TK ] J dQ W^T^Po-p, + p 3 + w) J , (5.6) 



r<£> (P ;ft+ft) 



9 2 1 



2iV c P° - 2m - 2Z r + 2^£ 

x 



{^V" ^ + X ) <5Q ' )(5r7 ~ 2A ^<*7<^T,] y StJ VF CTCT (w)W^ ff (ft +P4+LU) 

- [N c 5 aK S TX - S aX S TK ] J doj W aa (uj)T[ p ^ ri (P ;p 1 + p 4 + £J) J . (5.7) 

Fourier transforming to configuration space as in the previous section and with x, y representing the separations 
associated with ft + ft, ft + ft, respectively, the equations read 

rSUfl,;*) = (-O^ po _ 2w z 2 x r + 2 l£ {& c W + - 2N ^ 5 ^ w ^v? 

+ [N C S Sr.}w aa (y^[ c l rn (Po;y)}, (5.8) 

^(fl,;*) - (+^ po_ 2m _ 22r + 2 , £ [(^ 2 + - 2 T^(,) 2 

- [JVc^rA - 5 aA 5 T(t ] W^r^fb; x)} . (5.9) 
Further, we make the color decomposition: 

r&a = s aj s Tri r^ + s ari s TJ r^ 2 \ (5.10) 

where r( c ' p;1 ' 2 ) are now scalar dressing functions. Defining the functions 



\2N C J P Q -2m-21 r -i^ : {lTN c )W aa {x) + 2ie 
it is then straightforward to show that the solution to the above equations can be written 

2r( cl )(P ;y) = f+(P ;y) + f-(P ;y), (5.12a) 

2r( c2 )(P ;y) = f+(P ;y) - f-(P ;y), (5.12b) 

2T^(P ;x) = f-(P ;x) - f+(P ;x) = -2T^(P a ;x), (5.12c) 

2r(P 2 )(P ;x) = -f + (P Q ; X )- f_(P Q ; X ) = -2T^\P Q ;x). (5.12d) 

Putting the components together, we find 

T aT1TI (P ;x,y) = Si 7 S f Tr] [6 a7 5 T7l T^(P ; y) + 6 av S T1 T^ 2 \P ; y)] +6^6^ [s^S^T^ (I^; x) + S aTI 6 T1 T^ 2 \P ; x) 

= \ {(Sa-ySrr, + 5 arj 5 T1 ) [S } ai 8^f + (P ; y) - 8 s ari 8 s T1 f + (P Q ;xj\ 

+ 5 r ,)[5 f a7 5y.(Po;y) + S^f^x)] } . (5.13) 

As in the case of the qq systems, with this solution we return to the diagrams (f), (i) and (g). Writing out the 
explicit form of the energy integrals we notice that their form is identical to the quark- antiquark case, since the e 
prescription is similar, regardless of the internal quark (or antiquark) propagator. Thus, these diagrams are also 
vanishing. 
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Analyzing the pole structure of the above solution, Eq. (5.13), we notice that, as in the case of the qq systems, we 
have two different pole conditions (for either separation x or y): 



P res(x) -2m- 21 r - ± N c )W aa {x) = 



(5.14) 



which can be rewritten as 



Po res{x) =2m + g 



ctu D aiJ (uj) 



C F 



1±N C 

2N r . 



(5.15) 



With the plus sign (corresponding to the color antisymmetric term or /_ in Eq. (5.13) above), we find the condition 
for the finite pole position (i.e., finite integral even with infrared enhanced D aa ) P res (x), N c = 2, — 1; but with the 
minus sign (/ + , color symmetric), we obtain N c = —2, 1. Hence, the only physical solution, with N c — 2, corresponds 
to a color antisymmetric and flavor symmetric configuration, in agreement with our findings from Ref. [20] that 
bound states exist only for color antisymmetric SU{2) baryons (in that case flavor symmetry was implicit). Indeed, 
the physical poles in the flavor channels are the same (both correspond to /_), as demanded by the symmetry of the 
system. The homogeneous equation makes no reference to the flavor structure and can only have the flavor symmetric 
part. Also notice that whereas for the qq system there was a common unphysical pole in the normalization of I^ 4 ) or 
and then a distinct physical pole for the color singlet channel, in the diquark case each channel has a separate 

pole. 



VI. SUMMARY AND CONCLUSIONS 



In this paper, the Dyson-Schwinger equations for the 1PI and amputated 4-point Green's functions, for quark- 
antiquark and diquark systems, have been considered. At leading order in the heavy quark mass expansion and with 
the truncation to include only the (nonperturbative) gluon propagator and neglect the pure Yang-Mills vertices and 
higher order interactions, analytic solutions for the Green's functions have been obtained. We have found that the 
physical and unphysical poles disentangle, in contrast to the conventional phenomenological calculations. Indeed, 
this is one of the main results of this work, since in contemporary studies the interpretation of the spurious states 
appearing in the Bethe-Salpeter equation is still under debate. These results may hopefully be of further use in the 
construction of phenomenological models of mesons and baryons. 

The color singlet pole of the 4-point Green's function for qq systems leads directly to the linearly rising potential 
(confining energy for the bound state), for arbitrary number of colors, and in turn this can be related to the qq bound 
state energy arising from the homogeneous Bethe-Salpeter equation (considered in Ref. [20]). The second pole, which 
is situated at infinity as the (infrared) regularization is removed, is nonphysical and is part of the normalization; it 
is similar to the poles appearing in the gap equation and in the Faddeev vertex for baryons [21]. This simply means 
that the only relevant quantity is the bound state energy, i.e. the pole of the resonant component. In the case of the 
diquarks, the only physical solution, for N c = 2 colors, corresponding to a color antisymmetric and flavor symmetric 
configuration, again verifies our result from Ref. [20] that bound states only appear for SU(2) baryons and otherwise 
the system is not physically allowed. The second pole, corresponding to a flavor antisymmetric configuration, is not 
contained within the homogeneous Bethe-Salpeter equation (there, the flavor symmetry is implicit). We thus confirm 
that within the scheme studied here, the homogeneous Bethe-Salpeter equation furnishes exactly the right physical 
poles with no other solution and these are explicitly contained within the full Green's functions. 
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Appendix: Slavnov- Taylor identities 



In this Appendix we consider the Slavnov- Taylor identity for the quark-2 gluon vertex appearing in the text. The 
derivation is similar to the Slavnov- Taylor identity for the quark-gluon vertex presented in Ref. [20]. We start with 
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the full QCD action in the standard, second order formalism: 



S, 



QCD 



d xlq, 



tj°d + gT a j cr a +ij-V- gT a ~fA a 



m 



,1P 



4 p v 



where the (antisymmetric) field strength tensor is defined in terms of the gauge potential A® : 

F^ = d^-d v A; + gr bc AlAl. 
The action is invariant under a local SU(N C ) gauge transform characterized by the parameter 9 X : 

U x =cxp{-i9 a x T a } 

such that for infinitesimal 9 X , the fields transform as (with the notation a = A ) 

S aa = --d 9 a - f abc a b 9 c , SA a = -W9 a - f abc A b 9 c , 

9 9 
Sq a = -z8 a [T a ] a p q p , Sq a = i9 a q [3 [T a ] ^ a . 



(A.l) 

(A.2) 
(A.3) 



(A.4) 



The action Eq. (A.l) is invariant under gauge transformations and hence the functional integral Z is divergent by 
virtue of a zero mode. To surmount this problem we use the Faddeev-Popov technique and introduce a gauge-fixing 
term along with an associated ghost term. In Coulomb gauge (V-A a = 0), the new term in the action reads: 



Sfp 



d A x 



-\ a V-A a - c^V-D ab c b 



(A.5) 



where the gauge fixing condition is implemented with the help of the Lagrange multiplier A a , and c a and c b are the 
Grassmann-valued ghost fields. The spatial covariant derivative (in the adjoint representation) is given by: 



D ab = gab^ _ gf»cb A c 



(A.6) 



The action is invariant under a Gauss-BRST transform [1] whereby the infinitesimal spacetime-dependent parameter 
9% is factorized into two Grassmann-valued components: 9% = c%8\t- In Coulomb gauge, the gauge fixing term does 
not involve any time derivatives and hence one can define the infinitesimal variation 5Xt (not to be confused with the 
Lagrange multiplier A a ) to be time-dependent. The variations of the new fields read: 



Sc a = -\ a 5X t , Sc a = -\f abc c b c c 8X t , 5X a = 0. 
9 2 



The generating functional is given by 



with the source term 



S„ = 



Z[J] = J £>$cxp {iSqcd + iS FP + iS s } 7 



p a a a + J a -A a + cV + V a c a + + Q a X a + X a Qo 



(A.l) 



(A.8) 



(A.9) 



The Slavnov- Taylor identities in Coulomb gauge are derived from the observation that the Gauss-BRST transform 
can be regarded as a change of integration variables under which the generating functional is invariant. Provided that 
the Jacobian factor is trivial [8], we are left with an equation where only the source term varies: 



S [iSX t 



exp {iSqcd +iSfp + iS s + iSS s } 



<5A t =0 



- J 2?$ cxp { iSqcd + iS FP + iS s } J d*x5(t - xo) j-i (d° xP a x ) c% + f abc p a x a b c x 

-- g J? x V ix c a x + f abc J? x A b x c x - ix ax c a x T^q [jx - zc a x q^ aXax + - g X a xV a x + \f ahc rf x c b x ^ • (A.10) 
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The S(t — xq) constraint, appearing because of the time-dependent variation SXt, leads in principle to a non-trivial 
energy injection into the ghost lines of the Slavnov- Taylor identities (the ghost functions are, however, discarded in 
our truncation scheme). The above formula can be reexpressed by using the definitions Eq. (2.8) for the connected 
and proper functions, and after repeating the manipulations from Ref. [17, 20] one arrives at the following identity: 



= / #x5{t -x ) ji {dl <<>) c a x - r hc <<> [ 



' bc <<> \<ip h x irf x > +a b x c x \ 



V, 



(-v a 2 



<iA? > 



<icz> 



-f abc <iAl> tij (x) [<iJ h jx irf x > +A b x c x ] -A" <c> +\j abc <ic a x > [<iv b x vff x > +c b x c x ] 

9 z 



+*TZf, <iq ax > [<%Xp x vtfc> ~c x qp x ] +iT% a [<iXft x in> +c x qp x ] <iq ax >) ■ 



(A.ll) 



Note that functional derivatives involving the Lagrange multiplier result merely in a trivial identity such that the 
classical field A" can be set to zero [17]. Further, one functional derivative with respect to ic d must be taken, and 
it follows that the Slavnov- Taylor identities are functional derivatives of (as we will take further derivatives, all the 
fields and sources must be retained): 



= Jd 4 xS(t-x ) ^- l -{d° x <iai>)5{z-x)- f abc <ia a x > 

<iA1 x > <i? x ic d > -f abc <iA? x > U 3 (x) 



— d <ip b x irf x > -t<T b J dc S(z - x) 



— d <zJ* x ir}> -tA%S*°8(z-x) 



-iT% <iq ax > 



8id 



<iXp x vrf x > +S da 5(z - x 



j-t <iXf} x vrf x > S da 6{z - x)iq fj 



<iQ ax > 



(A.12) 



Starting with the above equation, the procedure is now to functionally differentiate with respect to the quark, 
antiquark and gluon fields. In principle, the resulting equation contains a large number of terms, however most of 
them simplify in the heavy mass limit and under truncation. To be specific, after taking the functional derivatives 
there are four categories of terms entering the equation — three of them are vanishing and one remains. Firstly, the 
terms multiplied by a spatial quark-gluon vertex do not contribute, since this vertex is suppressed at leading order in 
the mass expansion (as explained previously in Sec. III). Secondly, the terms containing a 4-point function Tq q Aa are 
also of C(l/m), due to the fact that in the corresponding Dyson- Schwinger equation at least one vertex in each loop 

term must be at tree- level (just as for the T qq A vertex), and we are at liberty to chose this to be (suppressed by 
the mass expansion) or a pure Yang- Mills contribution (explicitly truncated out). Thirdly, the ghost kernels arising 
from the functional derivatives also vanish, since they necessarily involve Yang-Mills vertices. Hence, under truncation 
and in the heavy mass limit, only the terms that involve a temporal quark-gluon vertex will survive. Explicitly, the 
equation Eq. (A.12), from which the Slavnov- Taylor identity for the quark-2 gluon vertex is derived, reduces to: 







d 4 xS(t - x )6(z -x)\— g (d° x <<>) + f aM <<> K - iT d fi <iq ax > iq [jx - iTf ja iqp x <iq ax > \ . (A.13) 



We now functionally differentiate with respect to iq ey , iq pt and and arrive at the following expression, after setting 
the sources to zero: 

= J dx a 5(t — x ) S(z - x) j--<9" <iqptiqeyio d x io e w > + f aed <iq pt iq Ey ta x > 5(x — w) 

+iT d e <iq p tiq ax i<J e w > S(x - y) - iT pa <iq ax iqeyZVw> $( x ~ *)} ■ ( A - 14 ) 

Identifying <iqiqioio> with T qqaa and Fourier transforming, one obtains (the momentum dependence of the right 
hand side is trivial) 



qqaa <l qqa ' qqa qi 

which under truncation (where T qqa — gives 



qqai 



(A.15) 



- qqaa 



jaedrpa _|_ ^X 16 ^ %T ± e T ] -|- j e ® ,ar £ a Q 



(A.16) 
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